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^ ■ Abstract 



A novel exact dynamical real space renormalization group for a Langevin 
equation derivable from a Euclidean Gaussian action is presented. It is demon- 
strated rigorously that an algebraic temporal law holds for the Green function 
on arbitrary structures of infinite extent. In the case of fractals it is shown 
on specific examples that two different fixed points are found at variance with 
periodic structures. Connection with growth dynamics of interfaces is also 

discussed. 
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I. INTRODUCTION 

Dynamical Renormalization Group (DRG) analyses of the critical dynamics of statistical 
mechanics models has a long history |I|^J^ with virtually all studies on regular translation- 
ally invariant structures. These approaches were, both in real and in momentum space, of a 
perturbative nature; in contrast, the principal aim of the present paper is the implementation 
of an exact approach, in real space, for the dynamics occurring on fractal structures. 

Consider a system described by a Hamiltonian (or action) if ({(yj}) with the field variables 
iPx defined on the lattice sites x. The simplest Langevin |l[] dynamics leading to equilibrium 
with the correct Boltzmann weight exp[— i/({(p})] is [^: 

at oLpx 

The temperature has been absorbed in the definition of H. The stochastic noise rix{t) is 
chosen from a Gaussian distribution 

V[{r,m) = N~' eM-Y. j dr^) (2) 

(A/"^^ is a normalization factor) which has a zero average and variance < r]r^{t)i]y{t') >= 
2D6,J{t-t'). 

We consider the simplest possible choice of the Hamiltonian H{{(f}), the Gaussian model: 

Hii^}) = J2 Y'fil - Y. 'P^'Py^ (3) 

X ■^ <xy> 

where a^ depends on temperature and the scale of ip^ is such that the second term in @, 
i.e. the sum over nearest neighbour sites, has coefficient 1. 

Besides being the starting point of more complicated models, the Gaussian model is 
related to many physical situation such as the properties of an ideal polymer in solution , 
diffusion processes and the growth dynamics of interfaces 0. Indeed if ax = k~^ for any 
X, then the hamiltonian (|^) describes equilibrium properties of an ideal linear polymer in 
solution with a monomer fugacity equal to k. In a^ = Zx, the coordination number of site 
X, the hamiltonian (^) can be rewritten as : 



H=lT.iV^- Vy? (4) 

and it is related to the diffusion process Icnown as the ant in the labyrinth 1^,0. 

If ifx represents the height of an interface above the substrate point x (in the sohd-on- 
sohd approximation), then eq. @ can be interpreted as the energy of an interface and the 
Langevin eq. (|l]) describes its dynamics. In a regular (hyper) cubic lattice this is known as 
the Edward- Wilkinson model for interface dynamics 0. 

We shall show, in specific examples, that in the presence of non-translational structures 
(such as fractals) even the simple model (|^) gives rise to interesting behaviours. Specifically 
we will show that two different fixed points are present unlike in structures with translational 
invariance where they collapse to the same fixed point. Part of the results presented here 
has already appeared in reference 0. 

The paper is organized as follows. In Sec.|I| we introduce the DRG technique for this 
problem first in the simple case of a one-dimensional lattice and then in a situation in 
which the couplings have a hierarchical structure. In Sec. ^ the non-trivial cases of fractal 
structures with uniform and non-uniform coordination number are analyzed. In the latter 
case it will be shown that there is an additional universality class. Sec. |V| contains both 
rigorous and heuristic arguments on a general network and for a non-linear case. The 
heuristic arguments are then checked by the numerical analysis of Sec. 0. Finally in Sec.|VI 
some closing conclusions are presented. 

II. ONE-DIMENSIONAL LATTICE 

We shall start with the simple case of a one- dimensional lattice. After a pedagogical 
example with uniform couplings, we will turn to a non-trivial example. In the latter case 
only DRG allows to obtain the asymptotic solution. 



A. Uniform couplings 

Let us start with a one-dimensional case in order to show how the method works. Equa- 
tion (P for the Hamiltonian (^ with a^ = a has then the form: 

^^ = (<^.-i(t) + <^.+i(t) - a(^,(t)) + r^^it) (5) 

It is immediately clear that if we were to choose uncorrelated noise, decimation would 
produce correlation between the nearest-neighbours noises. We thus assume a nearest- 
neighbour correlation to start with, that is: 

< r]^{ti)T]y{t2) >= 2D,J{ti - ta) (6) 

where D^^y = D^dx^y + Di5x±i^y As usual it is convenient to work in (time) Fourier space: 

a^^{uj) = (^^._i(u;) + <^x+i(^) + Vx{^) (7) 

where: 

|^e--/x.(^) (8) 

and a = a — iu. 

Upon decimation of (e.g.) the odd sites in favour of the even ones (i.e. solving eq. (^ 
for ip2x±i in terms of ip2x and <^2x±2), "we get an equation for the surviving (even) sites: 

/ 2 o^'h^M. ('^2x-2(a^) + 02x+2{uj)) , ; , . , . 

[a -2) = ^C2x{uj) (9j 

a a 

where we have defined a new noise Cx{^) as 

; / N ^ / ^ , (^2x-l(^)+^2x+l(^)) .T„x 

C2x(c^) = r]2xM H (10) 

a 

which is correlated as 

<4(^l) 4(^2) > = 2D^^y{27T)S{uJi + UJ2) (11) 

where as before D^^y = Do^^^y + Di6x±i^y and 



Do 




Di 


1 2 



(12a) 
(12b) 

Therefore the new noise keeps the same correlation as the original one and being a linear 
combination of Gaussian noises is itself a Gaussian noise. With the redefinitions: 



a' = a^ -2 



(13a) 
(13b) 
(13c) 



the Langevin equation preserves its original form. The static recursions relation for a is 
obtainable from ( |13b|) for a; = 0, i.e. 



/ 2 o 

a = a — 2 



(14) 



from which the recursions equation for u is derivable. Indeed if a = a — iu, then a' = a' — iuj' 
with: 



/ o -2 

uj = 2auj — tuj 



(15) 



The variance of the new noise is < ^z(i^i) Vy{^2) >= 2D' (27r)(5(to''i +^2) with D'iuo') 

. Do 
T{uj) ■ D{uj). Here we have defined the 1x2 matrix D 



^ 



and T is the transmission 



matrix which, in the long time limit u ^ becomes: 



Di 



T 



2a 
I2 



/ 



1 + 



2 4 



\ 



V 



2 

a J 



(16) 



It is important to note that the amplitude A, which is always allowed since it drops out 
from the final equations, is determined by the fixed point condition on D. The critical 
fixed point of equation ([M]) is a* = 2 (model (|^) in 1-d is meaningful only for a > 2). 
Since we are interested in critical dynamics, we will put a = a* = 2 in what follows. By 
setting the determinant of the system giving the fixed point for the matrix D equal to zero 



one obtains two solutions (two lines of fixed points corresponding to fixed ratios Dq/Di) 
associated respectively to A* = 2^/^ and A* = 2^1'^ . It is easy to see that the first choice 
leads to an unstable solution while the second results in a stable one. Indeed, let A+ and 
A_ be the two eigenvalues of the transmission matrix T and u+ and U- the corresponding 
eigenvectors. Then if we start with an intial state D : 

D = C+M+ + C-U- (17) 

(c-t being the coefficients of the expansion in the basis of eigenvectors) upon iteration of 
the RG transformation, the variance grows unbounded with the choice A* = 2^/^ while it is 
driven toward a stable fixed point D* = c+u^ with the choice A* = 2'^/^. 

We now have the means to compute the critical exponents of interest. From eq. (|13aj) 
we have, at the stable fixed point, for a generic scaling factor 6 > 1: 

0U^,D*)=b'f'0,{b'uj,D*) (18) 

where = 5/2, z = 2 is the dynamical exponent z and we have written 0{uj, D) instead 
of (^(2 — iuj, D). The exponent can be related to the scaling of the two-point correlation 
function: 

G^,y{tlM) = < ^x{tl) iPyih) > (19) 

where the average is over noise configurations. Using eq. (p!8D , we get: 

G,.,o(ti,t2)=62(^-^)G.,o(|,^) (20) 

The properties of the equal time correlation function can now be easily derived. If we define 
the function W'^{L,t) = Gl,o(^)^) we end up with the following scaling form: 

iy(L,t) = L"F(-^) (21) 

with a = (p — z = 1/2 and z = 2. 



Few comments are in order. If ip^ represents the height of an interface above the substrate 
site X, eq.(|^), with a = a* = 2 coincides with the growth equation proposed by Edwards 
and Wilkinson with unit surface tension z/, namely: 

^^ = z.VV.(t)+r/.(t) (22) 

where V^ is the discrete laplacian defined as 

^^^x = Y.i^y-^^) (23) 

Here notation y{x) means that the sum is restricted to the nearest neighbours of x. Then 
W{L,t) coincides with the width of the interface which, as a function of time and of the 
lateral extent of the substrate L, is found to exhibit the general scaling form W{L,t) oc 
L''f{t/L^), which reduces to W{L,t) oc t^ for t « L^ and to W{L,t) oc L" for t » L\ 
where a = /Sz and z is the standard dynamical exponent |1[]. 

For a (i-dimensional substrate Edwards and Wilkinson p found /? = max (^^,0j and 
z = 2 with an upper critical dimension of 2. Above d = 2 the interface is almost flat with 
a = j3 = 0. Our results on a one-dimensional case, thus coincide with the exact results. We 
shall now apply the same methodology to a non-trivial case. 

B. Hierarchical couplings 

A rather interesting case which lacks translational invariance and which can be solved 
only by RG techniques, is the case of a one-dimensional lattice where however the coupling 
between nearest-neighbours are hierarchically distributed as in Figure |T[ In the case of 
diffusion, those kind of models have already been studied |T^JTT||. 

In the present context we start from the modified Hamiltonian: 

HiW}) = Y. Y'PI - Y. ^x,yVxVy (24) 

where Wx,y = Wy^x have the hierarchical structure defined as follows (see Fig.|I|): 



Wx,x-i = Wx-i,x = { (25) 

where x = {2m + 1)2" (n, m = 0, 1, 2, ...) and Wx^y = if |x — y| > 1. The choice eo = 1 
can be made without loss of generahty by a proper rescahng of the time scale. 

As in the previous example we will work at criticality. This corresponds to the choice 

Oix 2^y=x±l '^x,yi 1-6. 

^ =n Y. ^x,y{¥^x - <^yY (26) 

^ <x,y> 

The implementation of the RG transformation closely follows the one used in the diffusion 
case [llO| , p!I[| and it is based on the decimation of the sites indicated by a circle in Fig.|l|. 
Interestingly, due to the particular decimation scheme chosen, the minimal set of parameters 
for the variance which are invariant under RG transformation, is D^y = D^dx^y + Di25y^x±i 
where one uses Di {D2) if w^^^ = cq (e„) respectively. In other words a single self-correlation 
but two different nearest-neighbour correlations are involved. The recursions turn out to be: 



(2 2\ 

^,x{uJ, D) = A^^^^^^,{u\ D') (27a) 



(«? - e?) 



e; = e„+i^^ i^ (n = l,2,....) (27b) 



uo' = ^ii±Mcu + 0(cu2) (27c) 



MU-') = ^.,(.) + ''"'"f,+ "■'"'"'"■ (27d) 

("1 -ei) 

where ai = —iu + 1 + ei and ULiyii) are left (right) sites which are decimated with respect 
to the barrier ei. Apart from the noise contribution, the recursions are the same as in the 
case of diffusion []TU|,|rT|. It proves convenient to define Rn = e„+i/e„ since, upon iteration 
of eq. (|27a|) , it flows to a fixed point R G (0, 1] and then the analysis of the fixed points 
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is reduced to a study of a two-dimensional map in the {ei,uj} plane. There are two fixed 
points 0= (e* = +00, cu* = 0) and A= (e^ = i?/(l - 2R),uj* = 0) whose stability depends 
upon the value of R. We will distinguish two cases: 

A) R < 1/2 : A is stable and O is unstable and y = dw = ln(2/i?)/ln2. The a; — ;> 0^ limit 

/ a{R) b{R) c{R)\ 



of the transition matrix T is T = (2/R)A 



c{R) d{R) e{R) 

V 1 y 



where a{R) 



2(1 - i? + R^), h{R) = 2(1 - R), c{R) = 2R{1 - R), d{R) = 2R, e{R) = 1-2R + 2R\ 
The value of A* which yields a stable fixed line is A* = (4/i?)^/^. Then, proceeding 
as in the previous example, one finds again a = ^^ = -^^ as expected since in this 
case df = d = 1. 

B) i? > 1/2 : O is stable and A is unstable and then y = dw = 2 as in the equal coupling 



case. The matrix T is T = AA 



/ ^ 1 -\ 



i 1 i 
22 



which yields A* = 2?!'^ (again as in 



V 1; 

the equal coupling case) and consequently a = 1/2. 



III. FRACTALS 

An interesting non-translationally invariant case, which can be analyzed with the tech- 
nique described above, is the case of fractal structures. We will consider several examples 
which are prototypes of different physical situations. 

A. Fractals with uniform coordination numbers 



The first example is the standard Sierpinski Gasket shown in Fig.^ whose fractal dimen- 
sion is df = In 3/ In 2. In this case, as on translationally invariant structures, there is only 
one fixed point. The renormalization group procedure is readily implemented in the Fourier 
transformed equation of motion (|lD for the hamiltonian (|]) which reads: 



ax'pxi'^) = Yl <^y(^) + Vx{^) (28) 

where a^ = ax—iuj and the sum is over nearest neighbours of x. The set up for the noise is the 
same as in the one- dimensional case. As before we put directly a = a* = A corresponding to 
the static fixed point. Then, following standard procedure, the RG transformation is carried 
out in two steps: 

a) The 0^ on the circled sites in Fig. ^ are eliminated in order to rescale the system by a scale 

factor 6 (= 2 in this case). This may be done by solving Eq. ( pSf ) for (fx{x = 1, 2, 3) in 
terms of 0y {y = A,B, C) (see Fig. g) and substituting the resulting equations back 
into (0). 

b) 0x and fjx are then suitably rescaled as follows: 

0bx{^,D) = A{uj + Q)0xW,D') (29a) 

Afi,,{uj') = {uj + 5){uj + 2)r),,(cu) + 2 ^ //^(c^) + (cu + 4) ^ fiy{u) (29b) 

y{x') yiix')) 

and 

u' = uj{5 — iuj) (29c) 

where primed (unprimed) sites refer to the surviving (decimated) sites under the RG 
transformation, and y{x) and y{{x)) are the decimated nearest-neighbours and the 
decimated next-nearest-neighbours respectively. 

The procedure for calculating the amplitude A* corresponding to the stable fixed point 
proceeds along exactly the same lines as the one- dimensional case: again the minimum set 
of parameters which are invariant under the RG transformation for the noise is {Dq,Di} 
corresponding to the self and nearest-neighbour correlation respectively. It is noteworthy 
that the presence of noise does not influence the renormalization of time, thus prompting 
the identification of the dynamical critical exponent z with the fractal dimension dw of the 

10 



walk on that structure giving the end-to-end distance R{t) of the ant in the labirinth at large 
t, i.e. R{t) ~ tV'^™, xhis is given by ( |29c| ) in the a; — ;> limit as oo' = b'^'^uj. One then finds 
the following critical exponents: 

, ^ 2 — (is 2 — ds , dw — df . . 

z = dw , P = — T — , a = — - — dw = 5 (30) 

with dw = In 5/ In 2 and ds = 2dp/dw = In 9/ In 5 is the spectral dimension describing the 
low frequency behaviour of the density of vibrational modes |]T^pr5| . A similar analysis can 
be carried out for a general d-dimensional Sierpinski Gasket. 

B. Fractals with non-uniform coordination number 

In the previous example we investigated the effect of the self-similarity of the structure 
on the scaling of the two-point correlation function for a Langevin equation induced by 
a Gaussian Hamiltonian. However in that example the coordination number was uniform 
(equal to 4) as in the translationally invariant lattice. The interplay between self-similarity of 
the substrate and non-uniformity of the coordination number has been recently investigated 



I^ . Contrary to what happens on periodic structures, it has been shown that on both 
statistical and deterministic fractal structures the equilibrium properties of (|^) are governed 
by two universality classes, corresponding to ideal polymers and to the ant in the labyrinth 



|T^ . We now study this effect on the Langevin equation on specific examples. The T-fractal 



(Fig.§) has Zx — 1,3 whereas the Branching-Koch-Curve (BKC) (fig. §) has Zx = 2, 3. 

We stress again that the appearance of two fixed points is a feature of non-uniform 
coordination number and a fractal structure. 

In the recursion relations we will derive below, the (minimum) set of parameters invariant 
under RG transformation is given by D^^y = Dq^ ii x = y {z is the coordination number), 
and Dx y = Di ii x and y are nearest neighbors and D^ y = otherwise. We will write 

I "A 
D = Dq3 and a = . With this distinction in mind, calculations along the 

\ "3/ 

same lines as before, yield the following recursions: 

11 



and 



V9ba;(a, D) = A Lp^{a , D ) 






ai= aia^ — 2 

3 - as/tti - 3 



(31a) 
(31b) 



a'n= o? 



(31c) 



In the above definitions f]y^^yo, represents the noise associated with sites of coordination 1 
and 3 respectively. 

In the static limit, cj = 0, one has a^ = cix and as found in |n|, the recursion equation 
( ^lc| ) admits two fixed points (see Fig|^): G = (a^, 03) = (1, 3), corresponding to the "ant in 
the labyrinth" , and P = (a^, Og) = (00, (1 + vT3)/2), describing the scaling behaviour of the 
ideal polymer [Q. (The symbols G and P denote growth and polymer respectively.) Indeed 
the linearized recursions (it is better to work with 03 and 0.3/0.1) show that G is unstable 
(two relevant eigendirections) with a thermal exponent ye = In 6/ In 2, which coincides with 



dw |T2||T^, the random walk dimension, while P is stable (one relevant and one irrelevant 
eigendirection) with yp = "^ ^^ — - related to the end-to-end distance R of an ideal polymer 
of length A^ according to A^ ~ R^^ 

Now let us turn to the dynamics. In the long time limit, i.e., u ^ 0, and near the above 
static fixed points one has a^ = o*. This implies that under renormalization u rescales as 
uj{b) ~ b^uj: thus the dynamical exponent z is equal to the thermal exponent y appropriate to 
the fixed point. From eq. ( |31b| ) one can easily calculate the new noise correlation functions 
near the fixed point. The recursion equation is: 



D' = T-D 



with T = QA- 



/ - - -\ 

' 2441 



3 5 1 

2 4 4 



/ 2o 



at the fixed point G and T = 2a A 



-1 A--2 



with a ^ 



\ 9 3 7/ 

\ 244/ 

i±^ at the fixed point C. 



(32) 



2a 1 a^ 
V 60 l + 3a/ 
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Thus, choosing the value of A so that the the largest eigenvalue is 1, one finds that the 



corresponding eigenvector D* is a stable fixed point. Near the fixed points from eq. ( pia| ) one 
readily deduces the scaling of 0x{^, D*) (we are writing ipxi!^, ^*) instead of (^^^(a* — zcj, D*)) 
in the u —>■ limit: 

ipU^,D*)=b''+'0x{b'uj,D*) . (33) 

From (|1|) and (H), one finds (6 = 2) 

ac = y^~ f = - at G , ap = ^^ ~ = 0.601.. at P (34) 

^ Zi Zi 

where c?/ = In 3/ In 2 is the fractal dimension of the structure. We will show below that 
the expression for ac is completely general. Since yc = dw and the spectral dimension, 
dsi describing the low frequency behaviour of the density of vibrational modes satisfies the 
scahng relation [0 ds = 2df/dw, the result (|34D for ac can be rewritten as ac = ^-^dw 



Taking into account that z = dw one gets /3 = (2 — ds)/4:. 

An extension of the T-fractals to sites with coordination z = 1,4 leads to exactly the 
same expressions as before with the same numerical value of ac = 1/2 with a slightly higher 
value ap = 0.678... for the other fixed point. Thus the ac is the same as in 1-d model of 
Sec. 0. 

A more interesting example because it is a structure with loops, is the case of the 
Branching-Koch-Curve (see Fig.|). Calculations in this case are rather tedious but quite 
analogous to the previous example and we will omit annoying details. Also in this case it is 
sufficient to start with a variance D^^y = -Doz^y.x + Di6\x~y\,i, where Dqz is defined as before 
and with two a's, ai and 02 corresponding to two types of coordinations. At the fixed point 
G=(a2 = 2,03 = 3) one finds again ac = {vg ~ df)/2 with i/g = dw = ln(40/3)/ln3 and 
df = In 5/ In 3 consistent with the previous claim of a general form. The other fixed point 
P= (og = +cx3,a3 = v5) yields yp = In 11/ In 3 and ap = In 2/ In 3 = 0.6309... Notice that 
in the above example ac < 1/2 whereas ap > 1/2. 
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IV. GENERALIZATIONS: GENERAL NETWORK AND NON-LINEARITY 

Eq. ([^) with Dxy = DS^y evolves to an equilibrium state with a probability distribution 
for {(p} given by: 

Veaiiv}) 0^ exp <^ — J2 ("f^ - Vyf \ (35) 



where x and y are nearest neighbor sites [|r3|. The equilibrium correlation function < 
{ipx — v^o)^ >T is defined as in eq. (^) calculated with the weight ( |53| ) can be shown |TT| to 



be proportional to the resistance VL^q between two fixed sites and x of the fractal network 
where conductances (equal to uj D) connect nearest neighbor sites that is: 

< {^^{t) - ipo{t)y >T ~ fi.,0 ~ kl^ (36) 

where the last equality defines the exponent ( for the resistivity. Scaling arguments predict 
that for fractals ( = dw — df (Einstein relation) [|16|. Since the roughness exponent a is 
defined by (see eq.(|19|,|0D < ((^^ - ^of >t~ |a;p", one finds a = ^^^ = ^^^1^. On the 
basis of the previous exact RG, we argue that the u renormalization (Eqs. (|3lD and (p3|)) 
is not influenced by the noise term in Eq. (|1]). We also note that Eq. (0) without the noise 
term is merely the diffusion equation with (p^ {t) interpreted as the probability of finding a 
diffusing particle at x at time t. If u renormalizes as lo' = b^uJ to leading order, then the 
mean square distance traveled after time t behaves asymptotically like t^/^ [|l^ implying 

Z = dy^ . 

We now show that the temporal behavior of the width of the interface W{t) for a substrate 
described by an arbitrary infinite network is also given by ^(2-^3)74 (^provided that ds can be 
defined 0]). This is a non-trivial result since we prove not only that the temporal growth 
has a power law form, but also obtain the exponent. The proof is heavily based on the 
rigorous results of Hattori et al. (HHW) [|1^]. We assume as initial condition ipx{t = 0) = 0. 
Let us go back to G^yit) =< ipx{t)ipy{t) > where the average is over the noise as defined in 
Sec.|. The formal (forward) solution of equation (^) can be written as: 
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roo 

Mt)=T. drU^yit - T)r]y{T) (37) 

y -"^ 

where Ux,y{t) is the retarded Green function, associated to (p2D , satisfying: 

(^ - vVl)Ux,y{t) = 6xJ{t) (38) 

that is: 

U.,yit) = 9it)ie'-''^\y (39) 

where 9{-) is the step function. 
Using (|38D , we have: 

G.,.,{t) = < ip,{t)ipx,{t) >=2DY, f dr a,,y{t - r)a,,,y{t - r) (40) 

y Jo 

roo 

= 2D [/.,.„(2(t-r)) 

•/ U 

where we have used the condition < rjx{ti)riy{t2) >= 2D5x,y5{ti — t2). 

Due to the initial condition Uxq,x{S^^) = Sxo,x, then Gx,xoit) is a solution of: 

i^^-ADu\/l)Gx,xoit)=5x,y2D (41) 

It is then easy to verify that a probability Px,xo (t) defined through 

Gx,xoit) = / drPx^xoir) (42) 

Jo 



is a solution of the diffusion equation 



[^^-uVl)Px,x,it)=0 (43) 



The results of Ref. |T8| can now be applied to Eq. ( |^ ) defined on the sites of an arbitrary 
graph. Since, under very general conditions, Pxo,xo{^) ~ t~'^°''^, from Eq.(|i2|) we get: 

W\L,t)\L=o. = Gxo,xo{t)--t^^ (44) 

where /3 = (2 — ds)/4 as expected from the previously solved cases. 
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From the same procedure one can recover also the exponent a. Indeed if we assume the 
standard ansatz which appears to be vahd for a generic fractal |n 



P.,.oit) = :^mx - xo\-n (45) 

and using eq.(|2|) and (|4|) one finds W{L, t) = L'^fit/L^) with a = dwi2 - ds)/4. as before 
(for ac and z = y). 

We now turn to non-linear growth on a fractal. We conjecture that the analog of the 
non-linear growth equation of reference PD[ is given by: 



+XY.Mt)-^^{t)f + r,^{t). (46) 

where now (fxif) is to be interpreted as the height of the substrate at position x at time 
t. The A term takes into account the lateral growth of the aggregate |^. When A = 
we recover Eq. (P^. Under rescaling of a factor h of the length, the left hand side and the 
first two terms of the right hand side scale as 6"~^, 6°~'^w' g^^d ^^a-dvy respectively, where 
]j-dw comes from the Laplacian on a fractal [^]. In the absence of the A term, dimensional 
analysis yields z = dw, (found above) whereas the presence of the A term leads to the time 
derivative on the left hand side balancing the A term: 

a + z = dw (47) 



with the u term being subdominant. Eq. (^) generalizes the exact equation p^ , p5| for 
Euclidean lattices which are characterized by d^ = 2. 



Additionally, following ref. p2|], if we identify the dimension of the noise term to be 
^-{z+df)/2 (^^]-^jg jg expected to hold in the linear case where the A term is absent) and 
extending the dimensional analysis, we find again eq. (|3^ ) for ac- 

In the presence of the A term, one may conjecture that the noise term scales as 5-(^+"'='c)/2 
where dc is the chemical dimension of the fractal [^. With this ad-hoc assumption one finds 



for equation (^61) : 
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leading to: a = 0.5064, (3 = 0.2789 in rf = 2, and a = 0.5170, /? = 0.2500 in rf = 3 Sierpinski 
gaskets. (For the gaskets in d dimensions dc = dj = \n{d + 1)/ In 2 and dw = ln{d + 3)/ In 2 

II)- 

We stress the fact that Eq. p7|) and ( ^8]) have been derived simply on the basis of 



power counting, and reduce to previously known approximations on translationally invariant 



structures 22,23 



V. NUMERICAL ANALYSIS 

In order to verify Eq. (|47| ) and check the conjecture leading to Eq. (^8|) , we have carried 
out computer simulations of growth models on Sierpinski gaskets. We considered a sequence 
of sizes of gaskets in both two and three dimensions, the largest one containing 1095 and 2050 
sites respectively. The value of a was estimated by comparing the scaling of the saturated 
roughness with system size, whereas /3 was deduced by studying the temporal dependence 
of the roughness for the larger sizes, (i.e. W{L,t) ~ L" for t » L^ and W{L,t) ~ t^ for 
t«L'). 

To check the correctness of the programs, we first studied the linear growth process in 
both d = 2 and 3 by adding and subtracting particles on the gasket sites with the same 
probability and found excellent agreement with the exact results of Eq. (|3^. 

The nonlinear case was then studied by carrying out simulations of the Kim-Kosterlitz 
||23| growth model for the gaskets; the numerical results confirmed the validity of Eq. (^7|). 
Specifically we found: a = 0.48 ±0.02 and (3 = 0.26 ±0.02 in ci = 2, whereas a = 0.48 ±0.02 
and f3 = 0.23 ± 0.02 in d = 3. 
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VI. DISCUSSIONS AND CONCLUSIONS 

Our discussion so far has been restricted to growth on fractal substrates with the particles 
arriving along a space dimension orthogonal to that in which the fractal resides. One might 
also consider growth on a self-affine substrate (with no overhangs, so that all sites are 
accessible to the incoming particles) with the growth occurring in the direction normal to 
the rough surface, but yet in the d-dimensional hyperplane in which the surface resides. In 
this case, the resulting growth would be characterized by the regular exponents in d — 1 
dimensions. This follows readily from two observations: 

i) The exponent a characterizing the growth is a measure of the equilibrium roughness 
and does not depend on whether the initial configuration is a self afiine or a planar interface. 

ii) It is reasonable that the exponent z which may be defined by studying the relaxation of 
small perturbations around equilibrium is typically the same that characterizes the approach 
to equilibrium from any initial configuration. 

The exponent /? is deduced directly from a and z. Thus the growth exponents in this 
case are trivially determined. 

Another observation is noteworthy. The results given here are valid for systems of contin- 
uous symmetry spins where, to the best of our knowledge, no analog of Henley's argument, 
valid for Ising spins undergoing Glauber dynamics, [^ is known. 



In conclusion we have presented a complete analysis of linear Langevin dynamics. The 
analysis has been carried out in real space and is based on a combinations of RG analysis, 
rigorous results and heuristic considerations. 

The primary results of our dynamical RG analysis are: 

1. For eq. (^21), i.e. for linear growth processes, the width of the interface grows 
algebraically with time and the exponent /3 = ^^^, where ds is the spectral dimension of 
the substrate. We stress that it is not common to explicitly demonstrate in a rigorous way 
such a power law behavior. Indeed, this result is valid not only for a fractal substrate, but 



for any generic substrate: ds can be rigorously defined [|T8| for almost [O] all infinite graphs 



consisting of a set of nodes and links joining sites that are defined to be nearest neighbors. 
2. For a fractal substrate, z = d]v, and a = ^^^^dw, where dw is the random walk 
dimension characterizing the asymptotic behavior of the root mean square distance, R, 
traveled after a time t, R r^ f^/d-w _ 
{dw, ds and the fractal dimension, df, are related by dg = 2df/dw WW^ )- 



3. For both fixed points of the hamiltonian @ the dynamical exponent z coincides with 
the thermal critical exponent y of the corresponding static problem. 

Result 1 is rigorous while 2 and 3 are based on exact arguments. For non-linear growth 
model (|46|) , heuristic arguments and numerical analysis suggests that eqs. (|48|j47|) are at 
least good approximations. 
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FIGURES 
FIG. 1. Part of a one-dimensional lattice with hierarchical couplings between the near- 
est-neighbours. Thicker bonds correspond to weaker coupling. Circled sites are eliminated after 
one RG step which scales the system by a factor of 2. 

FIG. 2. Part of an infinite Sierpinski Gasket in d = 2. Circled sites 1, 2, 3, ... are eliminated in 
favour of surviving sites A, B,C, .. after one RG step which scales the system by a factor 2. 

FIG. 3. Part of an infinite T-fractal in d = 2 which has coordination numbers z = 1,3. Circled 
sites are eliminated after one RG step and this scales the system by a factor 2. 

FIG. 4. Structure of the Branching-Koch-Curve in d = 2 which has ^ = 2, 3. Circled sites are 
eliminated after one RG step which scales the system by a factor of 3. 

FIG. 5. Flow diagram for the RG analisys for the T-fractal example. Fixed points G and P 
correspond to growth and polymer dynamic respectively. The line L is the critical line delimiting 
the unphysical region. 
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